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ABSTRACT 


This work deals with the correlation functions and their 
methods of obtention. A method without a "pure" time delay 
is considered which uses Laguerre functions type filters. 

The possibilities of utilization of the nonsymmetric Laguerre 
type filter are analyzed. A very inexpensive set of circuits 
which are useful to obtain the mentioned filters for low 
frequencies signals are presented, together with their 


computer analysis and laboratory realizations. 
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I. INTRODUCTION 


Correlation analysis has become the most powerful tool 
to separate the signal from the noise, especially in those 
applications where the signals are hidden beneath a blanket 
of other signals and extraneous noise. 

The development of instruments which can continuously 
sample even the most noisy signals and compute the correla- 
tion function simultaneously, have extended the use Of 
@errelation analysis into radio astronomy, fluid and solid 
state physics, neurology, seismology, geophysical exploration 
and other areas [Ref. 3]. 

Radar and sonar, of primary importance in military appli- 
cations, use correlation detectors because it has been shown, 
for the case of the binary decision problem, that no matter 
what the philosophy employed for the design of a receiver 
is, the best detector is always a crosscorrelator detector 
ef. 5]. 

Theoretically, an infinite signal to moise ratio can be 
mecained in the detection of periodic sigmais in noise 
through the correlation technique, but this would require 
an infinite time of observation. In practice, by using 
crosscorrelation practical measurements can be made in 
reasonable time. 

intctaliy> theoretical considerations @#@n correlation 


functions, their properties and methods to obtain them using 





pure time delays will be presented. Next, this thesis will 
present a method which does not use a pure time delay and 
yields results comparable to classical methods. 

The mathematical formulation of the method and the cir- 
cuits that can be used are studied and compared with the 
standard methods. Finally an approximation to the necessary 
type of circuits was analyzed by digital computer methods 


and checked in the laboratory using actual components. 





IIT. CORRELATION FUNCTIONS 


It is assumed that some physical process produces the 
time functions x, (t); x(t); Date y(t) simultaneously. 
It is assumed further that the physical process is stationary, 
in Other words it is not changing over time. Also it is 
assumed that the time functions are not zero and do not have 
ad.c. component or they have zero mean if the time functions 


correspond to a stochastic process [Ref. 3]. 


A. AUTOCORRELATION FUNCTION 

The time average autocorrelation function for determinis- 
tic or known signals is defined for two different classes of 
Signals, signals of finite average power (case II) and 
signals of finite energy or pulse signals (casel). 


For signals of case II the definition is: 


(22 (2) 2 lim 7 Py Cat exe (tetrad) dt (1) 
eo 1 1 


and for signals of case I: 


Q (2) 2 f x, (t) x, (t ae ej tele (2) 


=_—O 


If the time function under consideration is a stochastic 


process the autocorrelation function is defined as follows: 


R(t) = EL x, (t) x, (t + t)] (3) 





iemvnewstLOchasua.c process is ergodic 


Q(x) ment) (4) 


According to the definitions given above, Figure 1 shows 
an elementary device to obtain the autocorrelation function. 
The process implies multiplication of the signal by its 
delayed replica and averaging. If the instantaneous values 
of x, (t) and x, (t + t) are multiplied and the product. averaged 
Over a Ssuifieiently long time, the result is The autocorre- 
lation function if some means to vary tT have been added to 
Pew circuit. 

1. Properties of the Autocorrelation [Ref. 3] 

The autocorrelation function has the following 
properties: 
a. Maximum. 
The autocorrelation function is a maximum at 


t = 0, mathematically: 
IR¢ x) | <R(0) Schwartz Inequality (5) 


b. Value at Tt = 0 
ine value of the autocorrelatiom at 1 = 0 is 
connected to the total power of the signal. If x, (t) is a 
voltage, then for t = 0, ioe is the average power of the 
Signal (as measured on a 1 ohm resistor). Hence (k(0) is 
the average power of signals of case II ana(®(0) is the 


energy in the pulse signals of case l. 


10 





ss 





X(t) =COMERENT SIGNAL 4 
+WIDEBAND Noise 





RANDOM WAVE 2ERO AXIS CROSSING 
RATE K PER SEC., POISSON 
DISTRIBUTED 


x, (t) 





WIDEBAND NOISE 
RMS VALUE: Enq 
BANDWIOTH: & 


FIGURE 1: Typical Correlation Functions. 
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c. *unetdon 
The autocorrelation function is a function of 
Cheetorm of which is characteristic of the original 
mal x, (t). 
qd. Type of Function 





The value of the function for negative values of 
is identical to that for the same positive values. This 
ans that the autocorrelation function is an even function 


Ti 
Rox) AR (-1) (6) 


e. Averaging Time 
If the averaging time is long compared with the 
ciprocal of the lowest frequency in the original signal, 
‘a repeated measurements of the product for a given T will 
eld values very close to one another. 
me reriodic Functions 
isa x, (t) contains any periodic frequency compon- 
tse en +) will contain each of these frequency compon- 
Csmumeo ror certain time functions (R(t) is known if the 
iginal signal x, (t) is given. 
meee rc ysical Processes 
Any physical realizable process that produces 
time function like x, (t) is such that the value of the 
nection at time t becomes more independent of the value at 


+ tT as t gets larger. Thermal noise, and even the quantum 


2 





mechanical uncertainty principle, introduce randomness that 
cause a gradual loss of coherence between the signal and 
its delayed replica. This is true even for the oscillators 
in the most stable atomic clock. 

Hence for a signal arising from a real process 
the autocorrelation function approaches zero as t becomes 
sufficiently large. 

Dec ocme rence lame 

The value of t that causes a significant reduc- 
tion in the autocorrelation function is a measure of the 
coherence time of the original signal x, (t). 

i. The Frequency Domain 
A less obvious property of the autocorrelation 


function is indicated by: 


QQ (r)<34(w) (7) 


Equation (7) means that the autocorrelation 
function and the power spectral density are Fourier Transform 


pairs. Mathematically: 


o(w) 


QR (2) 


rRixe9*" at (8) 


x if o(w)eI#t dw (9) 


= 0° 


The power spectral density for deterministic 


ease JIL signals is: 


i 





o(w) = lim = KX pw) Ky mC w) (10) 


T'+-00 


NO 


- 


2 
Lim oq |X pl) | 


T'->00 
where the subindex T indicates that the function x, (t) age 
truncated and hence observed through a window of observation 
me duration T. 

The power spectral density for deterministic 


Base I signals is: 
% 2 
$(w) = X,(w) X, (w) = |X, (w) | (11) 


The same property holds for stochastic processes, 


in which case: 


R (1)<>6(w) ae) 
il o 
and ¢(w) = - sp EL |X p(w) |") Gig) 


j. White Noise 
The autocorrelation function for very wideband, 
uniform (white) noise with r.m.s. value En is an impulse 
function at t = O with amplitude EC. This means that one 
@aaracteristic of such kinds of noise is thet its instantaneous 
value is completely independent of the value at any other 
instant and that the coherence time of the process producing 


the noise is practically zero or very short. 
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B. CROSSCORRELATION FUNCTION 

While the autocorrelation function of a signal is equi- 
valent to the traditional technique of power density spectrum 
analysis, there is no classical analogy for crosscorrelation 
analysis. Crosscorrelation is concerned with the relation- 
ship between two different signals that arise in some common 
process [Ref. 3]. 

The mathematical expression of crosscorrelation is: 

T 


Caos = rae = a x, (t) x5 (t + ¢) dt (14) 


ieee GOpetcomon mune Crosscorrelation [Rer. 3] 


The properties of the crosscorrelation function are, 
in general, quite different from those of the autocorrelation 
Bane tLon. 

a. Type of Function 

The crosscorrelation function is not an even 


mumction Of Tt, which means: 


Cokes AR coe) | (15) 


However, KR on = 91 6T) (16) 
a relationship that has practical importance in obtaining 
Kon) for negative delays. 

b. Coherence 


imewerosscerrelavion function approaches zero 
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as t approaches infinity for signals that arise from real 
physical processes, due to the presence of noise and the 
mncertainty principle. 

Furthermore if x, (t) and x5(t) arise from two 


completely separate, unrelated processes 
Ro) = 0 ere) 


c. Frequency Domain 
The same as in the case of autocorrelation 


mumet Lons: 
(K (redo; 5(0) (18) 


In this case, however, the physical meaning of 
¢,5(w) is not so clear. It can be called the spectrum of 
etesscorre vation of the time functions x, (6) and x(t). 

| Figure 2 shows some typical examples of cross- 
Berrelation between functions. 

Crosscorrelation analysis provides a powerful 
analytical tool. The ability to measure the degree to which 
the signals that arise from a common physical phenomenon 
resemble each other as a function of the delay time between 
them, can provide a much deeper insight into the phenomenon 


being studied than a separate analysis of the properties 


of either signal alone. 
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FIGURE 2: Typical Crosscorrelation Functions. 
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IIIT. APPLICATIONS OF CORRELATION 


A. DETECTION [REF. 4] 
Given two functions consisting of a signal S(t) and 


noise N(t): 


X(t) = 8,(t) + N, (6) (19) 


X(t) = Sj(t) + No(t) | (20) 


it is possible to crosscorrelate them with the following 


result: 


( 


At 
yo(t) = lim ge oS [8 ,(t) + Ny (t) ILS, (t#r) + NQ(t+T)] at 


a Oe 


-R 


ip 
(21) 
eet) +K. Gs) +R, a") R, FD 
1-2 12 1-2 12 
All the terms of the above expression will be zero as 
16 approaches infinity. 
Now it is assumed X,(t) = X(t) and hence the four terms 


of the crosscorrelation become 
Rn R(t) Rot) Rn) #R ct) (22) 


In the expression thus eoeeneak is the autocorre- 


zi 
lation function of the signal itself and is a function of 
* other than zero even for large values of tT. The remaining 


verms approach zero as T becomes larger. So the first term 


18 


is a measure of the signal and the other ones are a measure - 
of noise. 

If now it is assumed that the frequency of the given sig- 
nal is known, an internal generator in the receiver can be 
used to generate the second input to the correlator and the 


two signals to crosscorrelate are: 


X, (t) S,(t) + N(t) (23) 


X(t) = Sp (t) (24) 


The above one is the typical case of a radar receiver, 


and if the two signals are crosscorrelated the result is: 


Ron Ro. +R ate) (25) 


S159 NSo 


Now the number of noise terms has been reduced to one 
and it is reduced to zero for large values of t. The pre- 
vious examples show how correlation techniques can be used 


to separate desired sigmals from the noise. 


B. ANALYSIS OF LINEAR SYSTEMS [REF. 3] 

It can be shown that it is possible to obtain a unit 
impulse response of a linear system by driving it with broad- 
band-white-noise and crosscorrelating this input with the 
system output, as is shown in Figure 3. The correlogram 
that results gives the same information about the system as 


if it were excited with an approximation to a delta function 
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and its output recorded on an oscilloscope or x-y recorder. 

This technique has important practical applications 
because of its immunity to internal system noise and because 
it permits to obtain the impulse response of a linear system 
in the presence of a signal in the linear system. 

With real time methods of correlation it is possible 
using the method described above to keep the response of 
critical systems under constant surveillance and to make 


optimizing adjustments. 


eu 





IV. METHODS TO OBTAIN CORRELATION FUNCTIONS 
There exist three methods to obtain correlation functions: 


A. METHOD THROUGH DEFINITION-TIME DOMAIN 
The easiest way to obtain a correlation seems to be the 
process indicated in its mathematical definition (1) and 


(14), this means: a delay t, is given to one of the signals, 


u 
then the other Signal and the delayed one are multiplied and 
finally the product is ononedpaine the result (1, ) one 
point of the correlation function. Varying t it is possible 
to obtain a graph representing the correlation function. 

A basic block diagram to perform the process was shown 
in Figures 1 and 2. The essential feature of such devices 
is the temporal delay tT between the two signals. When the 
Signals are fluctuating voltages or currents, the provision 
of the "pure" delay is rather a problem because the available 
delay devices produce distortion on the signals for long 
delay times. 

Another problem of this method is that it is essentially 
discrete, and the number of points of the correlation function 
that can be obtained in a given time interval depends upon 
the number of individual delay devices. 

Generally to avoid distortion, digital delay devices are 
used [Ref. 3] and hence it is possible to obtain a number 
@. discrete values of the correlation function without 


eastortion. 


2a 





Although this method is not cheap it is widely used for 


receivers which use correlator detectors. 


B. INDIRECT METHOD-FREQUENCY DOMAIN 

Equations (8), (9), (10), (11) and (13) set up the 
relations between correlation functions-time domain, and 
power spectral densities-frequency domain. 

It is seen from the mentioned equations that correlation 
of two signals, or one signal with a shifted version of 
itself, in the time domain is equivalent to a complex con- 
jJugate multiplication of their linear spectra in the frequency 
domain. 

In other words, the autospectrum and the autocorrelation 
are related through the Fourier Transform, the correlation 
being the inverse transform of the power spectrum. The same 
relation applies to the cross-spectrum and crosscorrelation. 

Although this method is the faster one it is expensive 
because to obtain correlations means have to be provided 
to obtain the power spectrums of the signals. Hence this 
method is typical of special signal processing machines 
which are designed to observe several characteristics of 
signals, and the power spectrums are not only used to obtain 


correlation functions. 


C. METHOD WITHOUT "PURE" DELAY-TIME DOMAIN 
This method has the advantage that a "pure" delay is 


not required. The correlation is obtained as a summation of 


e3 





terms and the method is based upon the properties of 
orthogonal functions [Ref. 2]. 

This is the only method essentially analogous and as 
will be shown in the following sections inexpensive 


correlators can be built for some special kind of signals. 
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¥V. CORRELATION WITHOUT A "PURE" TIME DELAY 


A. MATHEMATICAL FORMULATION 
Pr watalysis Problem 
It is supposed that any correlation mreseleral (Cay 
either autocorrelation or crosscorrelation, can be expanded 


in a series of orthogonal functions [Ref. 2], thus 


co 

(Ro) = Ea, 6 (t)Cw(r) I O< t<o (26) 
n=0 

where the 6 (t) are polynomials which form an othonormal set 

with respect to the weight function w(t) in the range 

teat <« ©. then the coefficients, a imeene series expansion 


can be found in the usual way and are given by 


n 


a= f Ror) e,(r)Lw(r) I at (27) 
0 


Now the system of Figure 4 is considered. The 
linear network has an impulse response h(t). So the filter 


output is: 


co 


x, (n3t) = af x, (t-u) h,(u) du (28) 


and the multiplier output is 


fee] 


x, (n;t) x, (t) = f x, (t-u) X5(t) h(a) du—s—« (297?) 


eS 





vaisaw 


SLIN313I4gaO0D JO NOIMNSZIGO -7 33anD\3 


NOMmoagNNoo Lo3aI0 
BAW I! LINW 


Oda 4inlgs 
sSv¥d MOT 
AYOMLAN 
(ziw)'x WVSANI1 


(3)°x 


(3)'X 
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being the output of the averaging device 


x, (nt) x(t) = ae x, (t-u) x5 (t) hu) du (30) 


pR(u) a (u) au (31) 
0 


Now comparing (31) to (27) it is seen that the 
coefficient an is identical with the average x, (n3t) X5(t) 


provided that the filter impulse response is so chosen that 
1-Y 
h(t) = 6 (t)Cw(t)] (32) 


The filters whoise impulse responses satisfy Equation 
(32) are called "orthogonal filters," and if they could be 
built, the coefficients an: would be obtaimed simply by 
switching in each of the filters in turn and determining the 
corresponding averages [Ref. 2]. 

Naturally, in any practical system only a finite 
number, say (N+1), of filters would be used. 


eeewovmonesis Problem 
Equation (26) can be written 


(R (1) 


Q(t) lw] (33) 


Rice a Ww a 


n=0 H 


or 


R (+) 


[w(t)]"* £ a h(t) Q<t<« (34) 


n=0 
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In particular, if a finite number of terms are 


taken, say N, it can be written as an approximation to (33) 


N 
RK ct) = [w(t)]°"? a h(t) > 


n=0 


If y is taken to be zero, (35) becomes 


N 
R ct) = Cw(t) > z a. h(t) (35a) 
n=0 
and if Y is taken to be 1/2 
(RQ N 
ny ft) = oe a, ae) (35b) 


Figure 5 shows the complete synthesis of the method. 
It is possible to use half the number of linear networks 
calculating first the coefficients due to the particular 
signals to be correlated, and then after introducing the 
values of a, the approximate correlation function can be 
seen on the screen of a CRT as a transient following the 
application of the input impulse [{Ref. 2]. 

A convenient system results if in fact the pulsing 
can be carried out repetitively at intervals that are long 
compared with the filter time constants, and yet sufficiently 
frequently to give a steady trace on the CRT screen. This 
can always be arranged in practice [Ref. 2]. 

It can be seen from Equation (35a) that if y is 


Baxen to be zero, the output of the circuit in Figure 5 has 
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to be multiplied by Putt] — to obtain the correlation 


function. 


B. TWO REALIZABLE SYSTEMS 

The applicability of all the preceding theory depends 
on realizing a set of filters whose impulse responses 
satisfy Equation (32). 

The polynomials 6 ft) and hence the impulse responses 
h(t) depend only on the weight function w(t) and the 
parameter jy. . 

The weight function w(t) = a4 e*t was chosen. 

1. Case of y = 0 

The transfer function H(s) of a simgle R-C cell 


(low pass filter) is: 





1 il fe! 
Shy SY) Sa (36) 
RC Toyi Sta 
iE 
being 
~i.lt 
RC 
Hac (s) corresponds in the time domain to the impulse 
response: 


hap(t) = a et (37) 


If now H(s) of a single C-R cell (high pass filter) 


is considered: 
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H Ss = + (38) 


cr‘s) = 


and the impulse response in the time domain is: 


h..(t) = -a e7% (39) 


cr‘ 
Now a number n of high pass filter cells are 
cascaded and between each cell a buffer amplifier is con- 
nected. Buffer amplifiers Pee aealiy infinite input 
impedance, zero output impedance and unity gain, and they 
isolate the cells. 
Hence the transfer function in the s domain for the 


system of n high pass filter type cells is: 


n 
H_ .p(s) = [H,p(s)]* = —=—— (40) 
nCR CR fe 3 ae 
rig a single low pass filter type cell is connected 
to the input of the circuit considered above via a buffer 
amplifier, the resultant circuit is of the form shown in 


Figure 6 and its transfer function in the s domain is: 





n 
= & a S 
heme uch “ste (|, > | CM} 
gl 
7 a s 
7 (s + q)otl sale 


Sil 





O=e t!1NDIID AJAL NOIMLSNNA BAYSNDOVI Ey SCIPS) [= 


oe te | 
‘ 4 re 





The inverse Laplace transform of the above one is: 


q® = 
a —n — e e oN (43) 
at ; 


which can be written as: 


ee 2 kT ka! 


yee L (at) 


at 


where ae’ is the weight function and L_ Cat) are the 


Laguerre polynomials and the product is defined as Laguerre 
fmametions . : 

The given expression for h(t) satisfies Equation 
(32) and the conditions of orthogonality [Ref. 2] and hence 
the circuit of Figure 6 can be used to build the correlator. 

2. Case of Y = 1/2 

Another filter based on the same weight function is 

the lattice structure shown in Figure 7. Im this case the 


s domain transfer function is given by 


ne) Fo: 
Wn) 
[ 
No] 2 
a 








H_(s) = ( ) (45) 


+ 
role 
we 
of 
of 2 


S) 


and the impulse response is: 


- & 
h(t) =% e e L (at) (46) 
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which is identical as the case of y = 0 except for a 
constant factor [Ref. 2]. 
3. Freguency Domain Considerations 

The nmonsymmetric filter, case of Y = 0, was used with 
good results to build an autocorrelator [Ref. 2] which used 
N = 10. Hence Ntl = 11 filters were used, being the first 
one a low pass filter (n = 0) and the rest combinations of 
alow pass filter and a different number of high pass filters 
isolated by cathode follower amplifiers. 

We saw in the ideal case, which:corresponds to 
perfect Jsolation buffer amplifiers, that the transfer 


function in the frequency domain was (41): 








» O S 
ence q ‘Ss aq: (47) 
and because of the isolation between cells it was taken as 
the product of the transfer function of a low pass filter 
and n R-C high pass filter cells. 

For the jw domain the above one becomes: 


oa 


JH#(ju)[° = 2% (48) 
(a 


A wey ott 
Meine a function of a and n. 

Figures 8, 9, 10 and 11 show the transfer functions 
in the frequency domain of a R-C low pass filter, an R-C 


high pass filter, ten high pass R-C cells cascaded and 
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ideally isolated, and of the ideal filters for n = 1 and 
n= 10, respectively. In all cases a = 103 was taken. 

The transfer funetions were obtained with the aid of 
the CORNAP circuit analysis program [Ref. 1] for the case of 
individual cells and for the case of multiple cells cascaded 
marame the product of the individual transfer functions. 

From the curves it can be easily seen which is the 
frequency range for the incoming signals and which the band- 
width of the system. Also if the incoming signals are given 
it is possible to see which mlteaae need postamplification. 

It has to be pointed out that the frequency domain 
considerations are only important in the determination of 
the coefficients (Figure 5) only case in which one of the 
input signals is convoluted with the impulse responses of 
mae fiiters. 

As the system is linear, by varying a it is possible 
to obtain filters for different ranges of frequencies. 

The correlator described in Ref. 2 used a = Q type 
filters to autocorrelate a sine wave which frequency was 
152 Hz, being the basic time-constant of the filter a = 843 
sec +. The results were practically identical to those 


which can be obtained with conventional methods. 


Ce. A SET OF PRACTICAL FELTERS 
The buffer amplifiers needed to build the filters have 
to have infinite input resistance, zero output resistence 


and unity gain and the system mentioned in Ref. 2 used 
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cathode followers as buffer amplifiers. A way to avoid the 
use of such amplifiers was studied in. this work. 
iMmemtacvor ©: Lsotation 

A filter of multiple cells can be constructed in 
such a way that the time constant is maintained the same 
through the cells and the resistance is multiplied by a 
factor of k from cell to cell as in Figure l2. It can be 
shown that this is equivalent to a certain amount of isola- 
tion between cells which depends on the value of the men- 


tioned factor: 


= —— = k = Factor of Isolation (49) 


The maximum value of k that can be used depends upon 
the necessary number of cells, and the values of the 
capacitors and resistors available. 

2. Computer Analysis 

In the Appendix a program is given which allows the 
calculation of the Laguerre functions for any values of a, 
nandt. Figures 13 through 23 show the Laguerre functions 
for a= 103 and for orders 0 up to 10. 

The CORNAP circuit analysis program was used to 
obtain the impulse responses of circuits of several number 
eeeecclis and different values of the factor of isolation. 

Impulse responses of circuits up to seven cells (n=6) 
were obtained because the results had considerable errors 


for circuits with a large number of cells. 
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The results showed that circuits with a factor of 
isolation of 10 have exactly the same impulse response that 
the ideal ones. The lowest factor of isolation that gave 
reasonable approximation to the ideal impulse responses was 
found to be 5. Hence, a factor of isolation of 10 can be 
used whenever it is possible and the use of values less than 
5 should be avoided. 

3. Laboratory Results 

Figures 24 and 25 show pictures of impulse responses 
Seecwe Kinds of circuits: eee for a = 103 and a factor of 
isolation of 10 for orders 0 up to 3, and second for the 
same a and k=5 for orders 0 up to 5. 

It can be seen that the practical results obtained 
agree with the theoretical considerations and computer 
predictions. 

It was not possible to build circuits with a larger 
mumber of cells because of the characteristics of the avail- 


able pulse generator, oscilloscope and components. 
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VI. CONCLUSION 


The realizability of the method of Laguerre type filters 
or without a "pure" time delay to obtain correlation func- 
tions was shown together with a set of very simple filters 
to achieve them for low frequency signals. 

The following things require further investigation: 

Multipliers and adder that could be used with the filters 
presented in this thesis to build a correlator as the one 
shown in Figure 5. | 

Factibility of building, through the techniques of Large 
Scale Integration, Laguerre type filters for several values 
of the base time constant alfa and with a behavior close to 
the ideal type of filter, which means that the circuits have 


to use buffer amplifiers. 
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APPENDIX 


PROGRAM TO COMPUTE THE VALUES 
OF LAGUERRE FUNCTIONS 


MAIN PROGRAM 


THIS PRCGRAM CALCULATES THE LAGUERRE FUNCT 
ORDERS O UP TON AS & FUNCTION CF THE ARGU 


THE PRODUCT ALFA*TIME IS TAKEN AS THE ARGU 
PASIC TIME CONSTANT :K ALFA 
NUMBER OF DISCRETE VALUES OF THE ARGUMENT :M 
INCREMENT OF THE ARGUMENT :SDELTA 

CIMENSICN eo 

READ 


= Wee 
m mo) 


AD(5,100)N 
READS) 100)M 
READ(5,1O1L)IDELTA 

STARTING VALUE OF ARGUMENT 
READ (5,101) X 
READ (5, 1LO2)KALFA 
L=N+1 
DO 3C I=1l,M 
CALL LAP(Y,X,N) 
DAMP=KALFA*EXP(-X) 
OC 20 J=1l,L 
Y¥( J) =Y(J)- DAMP 
WRITE(6,200) Y( J) 
CONTINUE 
X=X+CcLlIa 
CONT INUE 
FORMAT (12) 
FORMAT (F5.3) 
FORMAT( 17) 
FORMAT (F9.2) 
STOP 
END 
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